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Determinants whose elements involve random variables are discussed and expressions 
derived for the first and second moments. Applications are made to n-dimensional geometry, 
especially, to finding limiting probabilities for the event: "a, given point lies above a random 
hyperplane", under fairly general hypotheses. The random variable Anh—B is considered, 
where An and B are certain minors of the determinantal equation of the random hyperplane, 
and h \s a coordinate of the given point. 

An asymptotic expression for E{Al} is obtained, and it is shown that EIB^\ is of the 
order of (1/n) E\Al\. 



1. 

In this section real random variables X^j depend- 
ing on two indices i and j and having the following 
properties are considered: 

Let Ai, Aiy . . ., ^A he Z: random variables chosen 
from the A'^/s. Then we make the assumptions 
that 

(a) for any i a nd j, for any k, and for any choice 
of the AiS, the conditional mathematical expecta- 
tion 

E(X,j\A„ . . . , A,)=0, 

if Xij is not one of the Ai^. 

(b) For any i, j, for any k, and for any choice of 
the A/s, the conditional second order moment 

E{Xfj\A, . . . , A,)=<jI 

(where Ci is a constant that depends only on i) pro- 
vided Xij is not one of the AiS. 

Consider the following random determinant: 



A.= 



^21 



^12 

A.22 



. Xin 

' X2n 



X,a 



^nn 



We have the 



Lemma: E{A„) = andE{Al)^n\ n (^l 

i=l 

Prooj. Let Aij be the cofactor of X^j in A„; we 
can write 

E{An) = lZE{A,E{X,j\A,,)], 



From (a), E{Xij\Aij)^-^()] hence 

1 The preparation of this paper was sponsored (in part) by the OfTice of the Air 
Comptroller, USAF. 



Also 

A^ = IZAl^Xl + J:AjA,,X,jX,,, 

Now 

E{Y^AiiAijcXijXu] = 
JZE{A,jA,,X,jE{X,,\A,jA,,,X,j)} 

From (a) £'(A^u.|^i;,^ifc,A'i;) = 0; hence 

E[M]= j:e[ Ai XI } = S£;{ AiE{xi\A,;) : 

and from (b) this reduces to 

i 



(1) 



Now A\j is a A,,_i, which depends only on those Xi^ 
with ^5^L If w^e assume that the formula 



E{M]=n\j:a', 



(2) 



is true for 7i=l, 2, . . . , r, (1) shows that (2) is also 
valid for n = r-f-l; on the other hand (2) obviously 
holds for 71 = 1. Thus (2) is established by induc- 
tion, and the proof of the lemma is complete. 



dn= 



^11 ^12 

3^21 • 



'*'nl • 



^2n 



Nonrandom real numbers x^ and the determinant 
are now considered. 

Suppressing the k rows ii, 1^2, . . . , i^ in rf;,, a 
rectangular matrix is obtained with (n—k) rows; 
T[ii,i2, . . . , ik',(in] is called the sum of the squares 
of all determinants of {n — k)i\i order deduced from 
this matrix. 
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Let Dn be the following random determinant: 

a^ii + Xii a:i2 + Xi2 . . . Xm + Xm 

3^21"! ^21 • ... 



Z)„= 



Xn\\~Xn\ Xn2n~^n'. 



%~\-Xn 



Dn can be written as a sum of elementary determi- 
nants, the elements of which are either Xij or X^ but 
tiever Xij^Xij. One of these elementary determi- 
nants is dn] all the other elementary determinants 
are random determinants. Let S be one of these 
random determinants: S consists of k rows of 
XijS {k>l) and {n—k) rows of x^/s. For instance, 
suppose that the first row of S consists of the Xij] S 
is a linear homogeneous form in the Xi^'s, the 
coefficients of which are independent of the Xi/s. 
Thus according to (a) E(S) = 0; consequently 

Let ii, 12, . . . , ik he the k rows of S consisting 
of some Xi/s, that is, row I'l consists of the Xf^j, 
row ^2 of the Xi^j, etc. Employing the Laplace 
development of S in terms of these k rows ii, ^2, 
. . . , ijc, we can write S in the form: 



S=2*A,(iiJ2, . . . Jk)B(ji,J2, 



,jk), 



where Ajc{JhJ2, - - , jk) is a determinant of the 
preceding A type, the elements of which are those 
Xij with i € ii, 12, . . . , ik, j e ji, . . . , jk', 
BUuJ2, . • . J jk) is the algebraic complement of 
^k{ji,J2, • . . ,J7c); and the summation S* is extended 
over all combinations (ji,J2, • - , jk) of order k 
taken from the n integers 1,2, . . . , n. For two 
different combinations we have: 



E{Ak{JuJ2, 



, jk)-^k(ju 



Jl)]==0, (3) 



since there is at least one ja, ji say, which is not a 
j^, so that the product A^O'i, . . . , jk)'Ak(j[, • • • , 
jl) is a homogeneous linear form in Xi^j^, ^zgji? 
. . . , Xi^j^, and (3) follows from (a). Consequently 



E{S') = i:''B\j„ 
and bv the lemma, 



,jk)E[Al(ji, . . . ,jk)] 



E(S^) = j:%B%j,,.,.,jk)'k\u<rl) 



But obviously, 

i:*B\j„ . . . , j,)=T{H,i2, . 

Hence we have 



"^k'jdn)' 



,ik]dn)k\ n < 



Consider the product of two different elementary 
determinants Si and S2. There is at least one row 
consisting of some Xij that appear in Si and not in 
5^2, or in S2 and not in Si. For example, suppose the 
first row in Si, consisting of the Xi^'s, does not appear 
in iS'2; then the product S1S2 is a linear homogeneous 
form in the Xi/s, the coefficients of which are 
independent of the Xi^^s. Consequentl}^ by (a), 
E(SiS2)=0 and E{Dl) reduces to the sum of the 
squares of all the elementary determinants. This 
gives us the following : 

Theorem 1 ; Under assumptions (a) and (b) , 



E(Dn) = dn 



(4) 



E{Dl) = ^kl(^^niiM, . . . ,ik;d^) n <). (5) 

In (5) the summation S* is extended over all combi- 
nations {i\, 12, . . . ,%) of order k of the integers 1, 
2, . . . ,n] for convenience we put 



r(i,2, 



,n]dr) = l. 



An interesting feature of formulas (4) and (5) is 
that these formulas do not depend on the probability 
laws of the X^/s. This fact remains valid even if the 
(t/s depend on j, but (5) becomes more complicated; 
on the other hand it is possible to compute the higher 
moments of D^ by using similar reasoning and suit- 
able assumptions. However, the formulas seem quite 
complicated. It also appears difficult to obtain the 
probability law of Dn, even under such hypotheses as 
that the Xij are normally distributed. 

Of course if 0-1 = 0-2= . . . =(Tn=^<^, where c is a 
constant independent of i, (5) reduces to : 



where 



E(Dl)=J^ Tk(dn)k\a'\ 

k = 

Tk(dn) = ^*nii, . . . ,ik;dn). 



(5a) 



From (5) certain interesting results may be deduced, 
which are perhaps already known. These results are 
connected with the geome trical interpretation of the co 
efficients of (5). Let £'^ be an 7i-dimensional Euclidean 
space, with orthogonal coordinates. The n numbers 
Xii, Xi2, . . . , Xin may be regarded as the coordinates 
of a point P^ in En, and Xn-^-Xn, Xi2-\-Xi2, . . . , 
Xin-\-Xin can be considered as the coordinates of a 
random point Mn in E^. That dn is invariant under 
orthogonal transformations is well known. If 
(0,PiP2 . . . Pn) is the volume in E^ of the paral- 
lelepiped formed by the vectors OPi, OP2, = . . , 
OPn we have : 



dn=±{0,PiP2 



Pn)^ 



An analogous interpretation holds for Dn and 
(0, M1M2 . . . Mn); hence E{Dl) is invariant under 
an orthogonal transformation. If the XijS are nor- 
mally distributed, this property also persists. The 
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(Ti remain unaltered also. We conclude that the 
coefficients T{ii, . . . , 4; ^J are invariant. 

T{ii, . . . ,4; <^ J depends only on P,-^, . . . ,P,-^; 

consider the set of points 7^ p T^i 2' • • • ? P^i^ defined 
as follows : 



op;,=op,^+x2iOP,, 



op; =op. +x,iOP. + . . . +x,.,_iOP, 



k-V 



where the \ij are such that the vectors OP/^, O/^/^j 

. . . , OP;^ are mutually orthogonal. From a classical 
property of determinants we can replace P^^, . . . , 
^h ^y P'iv ' ' y ^h without altering T(ii, . . . ,'i^; dn). 
The invariance of T{ii, . . . , ifc] dn) under orthogonal 

transformations implies that the vectors OPi^(a=l, 
. . , , k) can be taken as coordinate axes. Then it 
becomes obvious that 

r(ii, ...,{,; ^,):=(o,p;^...p;/, 

the volume (0, P-^ . . . P-) being considered as a 
volume in a /c-diniensional subspace of E^; but it is 
also obvious that 



(o,p;...p;/=(o,P..P., 

and we have established: 
Theorem 2: 



■ p^J 



T{i„... ,H;d„)={0,I\. 
and consequently, 



p^.y 



(6) 



(7) 



On the other hand, considering Pi,. . ., P„ as fixed 
points in En, and 1] as a moving point, (Q, P,- ... 
P,^y and 2* (12, P,^ . . . P,^y are given by (6) or 
(7) under a simple change of the origin of coordinates. 
If one puts 



and 



(n,P,^. . .P,X=C 



(8) 
(9) 



where C and C* are any positive constants, then eq 
(8)and(9)definequadricsQ(i^,. . ., i^; C) and Q^{C*), 
respectively, of elHpsoidal type. Thus it follows 
from (5) that the relation 



E{{^,M,. . .MnY}=e, 



(10) 



where e is some positive constant, also defines a 
quadric Q {e) (which is in general an elHpsoid). 

2. 

In Euclidean ?i-space with the coordinates defined 
as above consider n fixed points Pi, Po,. . ., P^, and 
a fixed point S on the axis Oxn with coordinates 
(0, 0,. . ., 0, h). Let Ml, M2,. . ., Mn be n random 
points, and consider the random plane H, (i, e., an 
(?i-l) -dimensional linear variety) defined by (Mi, 
M2,. . ., Mn). The equation of n has the following 
determinantal form 



U= 



Xi X2 

Xii-|-Aii Xi2-\rXi2 . 



^1, «-l~h^l, w-l -^In 1 



^nl~T^n 



-l~r-^n,w-l -^nn 1 



Developing this determinant in terms of the elements 
of the first row and their cofactors, one obtains 

U=f^AjXj-B=0. (11) 

For convenience we set 

Yij=Xij+Xij (i=l,2,. . .,n;j=l,2,. . .,n—l). 

Let Mi be the point (F,:2v • • Yi,n-i, 0) and V the 
linear variety of ri — 2 dimensions defined by 



= 0, 



il=0, 



(12) 



and let Hi be the orthogonal projection of Mi or 
Mi on V. An elementary computation shows that, 
in absolute value, 



M:Mi__ \Xi 



H.M'i HiM[ 



Vn-l 
\An\ 



(13) 



We now suppose that: 

(A) the distance between P^ and P^- is equal to 1, 
for any i and j. In this case the n points Pi form a 
regular polyhedron in the variety ^^^=0; this poly- 
hedron admits l^\ linear varieties Vij as varieties of 

symmetry, Vij being defined as that (ti— 2) -dimen- 
sional variety formed by the points in Xn^O which 
are equidistant from P,- and Pj (ir^j). Now suppose 
that F is a moving (n— 2) -dimensional variety, and 
consider the quantity 

s=max (distance from F to P^). 
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s attains its minimum value when Fis a Vu, in which 

case s=—' Consequently in all cases there is at least 

one point Pf such that the distance from V to P^ is at 

least -• 

Now suppose that n tends to + oo with condition 
(A) always holding; we make the following assump- 
tions:' 

(B) lim pAm^x \X,n\<e\-=l for any €>0; 

(C) there exists a positive number a<^— and inde- 
pendent of n such that: 

limPrJMaxi]^,<a2 

n-*co ( i ;=1 



:1. 



If for instance the distance from P^ to y is a maxi- 
mum when i = lj then except in cases of small 
probability: 



HiM[>--a, 

and 

Thus we have: 

Theorem 3; Under hypotheses (A), (B), and (C), 



lim Pr I 



iir>'') 



-0. 



Remark: Assumptions (B) and (C) are not particu- 
larly restrictive. They are satisfied, for instance if 
the Xi/s are independent and normally distributed 
with 

E{Xf,)^^, (14) 

p being<- and independent of i, j, n; also when the 

Mi are independently distributed with uniform 
probability density over the interior of the sphere 

of center Pi and radius p<C^* 

Denote the center of gravity of Pi, . . ., P^, by G^, 
with coordinates (fi, . . ., f„) given by 



n i=i 



(i=l,2, . . . ,n), 



(15) 



^„ = 0. Let K be the center of gravity of Mi, . . 
M„ with coordinates ^j+Zj (j=l,2, . . ., n). Put 



1 ^ 

n i=i 



(16) 



We make the following assumptions: 

(D) the distance OG remains bounded as 7i->+ °° ; 



(E) the Xij satisfy assumptions (a) and (b) of 
section 1, with 

-y/n 

where p^ is a constant independent of i and n. The 
algebraic distance L from G to n is equal to: 



(17) 



71-1 

K belongs to IT, so that 

\L\<GK. 



But the coordinates of the vector GK are the Z/s 
defined by (16). From (E) it then follows that 

E{GK')<^' 



Hence GK and a fortiori L tend toward zero in prob- 
ability. Write (17) as 





B A„ 


Vs^' 


An 1 n-l 

n-l 

1 •'■=1 



Vn-1 
' + TAJ 



We have by the Schwartz inequalit}^ 

n—l I In— I In — ] 



SO that 



Now 





V 1 / 



zero in probabihty. Therefore under the assump- 
tions (A)— (E) we have: 

Lemma: Asn^<^, \B/An\ tends to zero in prob- 
ability. 

The algebraic distance fromS to IT is 
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hA„-B 



hAn-B 



\ \An\ \An\J 



B 



,v: 



ii-l 

'JAJ 




(18) 



Suppose that A->0; then S is above n if AA„>0,or 
what is the same thing, {hAl—BA„)y-0. We may 
write this last expression as 

Thus according to the above lemma (and under the 
same hypotheses, viz., assumption (A)-(E)): 

Theorem 4. IJ /^>0, the probability that S is above U 
tends to 1 when n->+ 00 ; ifh<iO, this probability tends 
to zero. 

The case //=0 is not clear-cut. But if, for in- 
stance, the M/s are independent with the same con- 
tinuous probability law which is symmetric about 
the plane Xn = ^, then clearl}^ Fr { S > f } == | . 



Keeping the notations of section 2, we assume that 
P^Pj=l for any i, j and that the Xt/s satisfy con- 
ditions (a) and (b) of section 1 with 0-1= .. . =(t^ = 
p/V^ where p is a constant independent of i and n. 
The other assumptions, (B), (C), and (D) of section 
2, will not be required. We wish to obtahi some in- 
formation about the random variable: 



H=Anh-B, 
Theorem 1 implies immediately that 

E{H)=VnK 

where Vn is the determinant 



(19) 
(20) 



V„-- 



X21 



^' n\ 



X\2 ... 3^1, n-\ 

X%2 ' • • ^2, n-1 



Xn2 • • • -^n, n—1 

Vn, as w^e have seen, is given by 

That is to say (cf. formuhi (33) of the Appendix) 



F»=±-y'2^^i* 



(21) 



Let w„ be the determinant 

^12 



^11 

a: 21 



• ^2,w-l 



Xnl • • • Xn,n —1 O 

For E{W) we have at once: 

E(H')=h'E{Al)+E{B'). 



(22) 



Since B is a determinant of the form Dn (section 1) 
formula (5a) gives the result that 



E(B')=j:T,iw„)k\^ 



k=l 



(23) 



An is not exactly of type Dn, but in entirely similar 
fashion we find 



E(Al)=^Tl{Vn)k\^- 



(24) 



Tl(Vn) is the sum of the minors of order {n — k) of 
Vn obtained without suppressing the last column of 
Vn- For instance 

n{i,,...,w,v.)=(Pj„p^,Pj,...Pj,^_,y, 

where ji,J2,. . ',jn-k run through all the numbers 
1, 2, . . .,71 except 'ii, 12, . . • , in-k- Hence 



and 



J- k\^l, • • • j'^Aj ^n) — y n-k — 2n-k-2 



)-n+A;+l 



Substituting this expression into (24) one finds 



(25) 
(26) 



^2^S n^ ^^^^• 



(27) 



If we suppose that 2p^<l, (27) yields immediately: 



E{Aiy 



n 



'l_2p2-2«-i' 



(27a) 



as n-^ 00 . 



B is of the form D„; we know that ^^(C*) is a 
quadric, and it must have as planes of symmetry the 



(^j planes of symmetry Vtj of the polyhedron 
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(Pi, P2, . . ., Pn)> Consequently the center of Q,(C^*) 
is G, and QkiC"^) is a sphere since I j > n, and 



we can write: 



T,{Wn) = Un)OG' + mWrd, 



(28) 



where Tt{wn) is the value of Tj,{Wy^ when OG=0. 
On the other hand 

T,(ii, . . .,ijc;wn)=(0,Pj^Pj^ . . . P,_;, 

where 7*1,^2, • • • , jn-k are all the numbers 1, 2, ... , 
n except ii, ^2, . . . , ik- If 0= 6^, we find (cf . appendix, 
formula (36)): 



(1/n!) \Vn\ of Pn in En-i- Let G^i be the center of gravity of 
the points Pi, P2, . . . , Pi', hn the distance between (zn-i and 
Pn', and dn the distance between G„ and any P^. We have 



\Vn\=hn\V„ 



(n-1) GnG^.-l + GnPn = 0. 



(31) 



Gn-h Gn, Pn are on the same ''straight hne'' (one-dimensional 
linear variety); and PiGn-iPn is a right triangle; hence: 



It follows that 



1 = PiPn - P1GI-, + (?._iP; 






r*(ii, . . . , ij] Wn)^ 



?i2" 



and 



Hence 



n(».)-(j)i^ 



;JT(».)«?^=;it t '"-')-J-'-+') (V)-. 



'/i^ n2^ fr{ 



If 2p^<l, (29) becomes asymptotically 



12TKwn)kl 



1 



Ip^ 



yt=i 



n' 7^2^ (l-2p')' 



(29) 



(29a) 



We can now compute the value Tt'^iWn) of Tk(Wn) 
for Pn=0 (cf. appendix, formula (38)). Writing 

Tr(Wn) = Un)GPl + mWn), 
we have an equation which determines X^^ (n) ; we find 



and 



. ._ n(n-2)(n-3) . . . jn-k) 
M^j- (A-l)!2'^-*-i 



j:Un)k\- 



1 ^ n(n-2)(n-2) . . . (n-k) 

(30) 



hl=i- 



(J^D" 



/»„^-i. 



(32) 



Clearly h^ is known; the recurrence equation (32) can be 
solved for hn: 



L l-JL 



It follows by recurrence from (31), IF3I being equal to V3/2. 
that: 



On the other hand, 



and 



Y 2»-' 
1 /re^ 1 

777.2 _ 1 _!/ " i\ 



(33) 
(34) 
(35) 



We now compute the volume {Gn, P1P2 . • • Pn-k)- If Lnk 
is the orthogonal projection of Gn on the linear variety 
(P1P2 . . . Pn-k), we have: 

{Gn,PlP2 . . . Pn-kr = Gj:ik\Vn-k\'. 

But Lnk = Gn-k, so that 

GnLnk^^GnGn-l'TGri-lGn-2'T- • • • ~r Gn-k+lGn -k' 

As in (35), 



and so 






■/c)' 



/c = l 



Asymptotically this expression becomes, for 2p^<l, 

n ^2k I 2p2 



j:\k{n)k\ 



v} 2^-i(l-2p7 



(30a) 



We verify that E{B'') is of the order of {lln)E{Al), 
which agrees with the lemma of section 2. In the 
particular case when G=0, E(B^) is of the order of 
(l/n^EiAl). 

Appendix 

Given n points Pi, . . . , Pn i n an (n— I) -dimensional 
Euclidean space En-i such that PiPj=l for any i and j 
(^Vi) (that is, these points are the vertices of a regular poly- 
hedron Pn) we wish to find an expression for the volume 



((?n,Pl . . . Pn-k)'=- 



k 



n2"-* 
We can also compute (P„, P1P2 • . . Pn-k)'- 

iPn,PlP2 . . . Pn-k)' = PJ^lk\^--^\' 



(36) 



PnLnk — PnGn-l~\~Gn-lGn-2~^ 

4('+^)' 



-Gn-k+lGn- 



Therefore 



{Pn, Pi . . . Pn-k)' = 



(n-/c4-l) 



2n-k 



From (37) it follows that 

Los Angeles, January 11, 1951. 



(37) 



(38) 



470 



